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Abstract 

The value of the safety factor on the magnetic axis of a finite-beta spheromak is shown to be 
a function of beta in contrast to what was used in P. M. Bellan, Phys. Plasmas 9, 3050 (2002); 
this dependence on beta substantially reduces the gradient of the safety factor compared to the 
previous calculation. The method for generating finite-beta spheromak equilibria is extended to 
generate equilibria describing toroidal magnetic “bubbles” where the hydrodynamic pressure on 
the magnetic axis is less than on the toroid surface. This ” anti-confinement” configuration can be 
considered an equilibrium with an inverted beta profile and is relevant to interplanetary magnetic 
clouds as these clouds have lower hydrodynamic pressure in their interior than on their surface. 
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I. INTRODUCTION 


f. o, 


In Ref. [1], one of the authors (PMB) examined analytic forms of finite f3 spheromak 
equilibria and used a well-known expression for the value of the safety factor q on the 
magnetic axis, denoted as q ax i S , to argue that finite f3 causes the beneficial effect of a much 
larger q gradient than when /3 — 0. However, co-author (RP) numerically calculated q ax is 
for these finite (3 analytic equilibria and found numerical results substantially different from 
the q ax i S given in Ref. [1]. The reasons for this difference are identified as resulting from a 
subtle misuse of an expression for q ax i S . Resolution of this issue revealed that the analytic 


equilibria presented in Ref.[l] could be extended to give an interesting toroidal equilibria 
where the pressure on the magnetic axis of a toroid is lower than the pressure at the surface 
(edge) of the toroid rather than higher as in a tokamak; i.e., the beta profile is inverted and 
the configuration is bubble-like. Increase of a parameter 7 (to be defined below) results 
in solutions to a Grad-Shafranov equation evolving from characterizing finite f3 spheromak 
equilibria, to a conventional zero /3 spheromak, to magnetic “bubbles” which are tokamak- 
like configurations having inverted /3 profiles, and then to a tokamak with conventional f3 
profile. This evolution is characterized by the ratio of two Bessel functions changing sign as 
their argument 7 is progressively increased. Interplanetary magnetic clouds are an example 
of the magnetic bubble situation because on the magnetic axis these clouds have lower 
hydrodynamic pressure than at their edge. Magnetic clouds have been previously modeled 
using numerical solutions to Grad-Shafranov equations [ 2 ], [ 3 ] in a slab approximation (i.e., 
equations are solved in Cartesian geometry in the xy plane with the z direction ignorable); 
the model presented here differs by being analytic and axisymmetric (i.e., equations are 
solved in cylindrical geometry in the rz plane with the </> direction ignorable) so that, in 
contrast to a slab approximation, toroidal geometry effects are inherently included. The 
analytic model has only a few parameters and so has less freedom than a numerical model 
but nevertheless has the useful feature of revealing parametric dependence and scaling. 
The analytic model also offers the possibility of providing a useful framework for other 
calculations, for example, calculating particle orbits in an axisymmetric cloud; the virtues 
of developing a repertoire of analytic solutions to the Grad-Shafranov equation has been 


discussed in Ref. 


4 ]. 
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II. BASIC RELATIONS 


We use a cylindrical coordinate system {r, 0, z} and consider the general axisymmetric 
magnetic field 

B =-^-V0 x V0 + Tb,rV0 (1) 

Z7T 

where 0 is the poloidal flux function and B$ is the toroidal field. The 0 direction is called 
the toroidal direction and any direction lying in the poloidal plane (rz plane) is called a 
poloidal direction. From Ampere’s law the associated current density is 


We are interested in configurations where the poloidal flux function has a local extremum 
in the r, z plane; both spheromaks and tokamaks are this type of configuration. The location 
of this extremum is called the magnetic axis and its vertical location defines the z origin 
while its radial location is defined to be r ax i S ] 0 is thus at a maximum or a minimum at 
r — Taxis, z = 0. If 0 is at a maximum on the magnetic axis then J$ is positive at the axis 
whereas if 0 is at a minimum on the magnetic axis then J^ is negative at the axis. 

Spheromaks and tokamaks are characterized by the safety factor q which is the number of 
times a field line goes around toroidally for each time it goes poloidally around the magnetic 
axis. Tokamaks typically have near-unity q on the magnetic axis with q increasing with 
increasing distance from the magnetic axis whereas spheromaks have near-unity q on the 
magnetic axis and q decreasing on moving away from the magnetic axis. The gradient of 
q, denoted as q ', provides stability properties and detailed calculations show that a zero 0 
spheromak has small q'. 

The safety factor at the magnetic axis is given by jsj] 


Qaxis 


T/2 + e -l/2 £ 


<j),axis 


/^0 


( 3 ) 


where 



is a measure of the cllipticity of 0(r, z) in the vicinity of the magnetic axis such that e > 1 


corresponds to vertically elongated equilibria (prolate) while e < 1 corresponds to vertically 


shortened equilibria (oblate). The force-free relation /i 0 J^axis/B^axig = A was invoked in 
Ref. [51] to give q ax i S = (e 1//2 + e^ 1 / 2 ) /(A r ax i S ) but this result is valid only if the plasma is 
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indeed force-free (i.e., has zero f3 and equilibrium given by V x B = AB). If (3 is finite, then 
Po J 4 ,,axis 7 ^ A B^ axis and it is necessary to calculate the actual value of no J<t>,axis/ B^ axis by 
consideration of the details of the finite (3 equilibrium. 

To do this, we start by defining /3 re i 


drel hO 


Pn.xi<! -Pi 


Ic 


B, 


( 5 ) 


where P ax is and Pi c are respectively the hydrodynamic pressures on the magnetic axis and 
on the last closed flux surface. Positive /3 re i thus corresponds to a conventional (3 profile 
whereas negative f3 re i corresponds to an inverted (3 profile. This definition differs from that 
used in Ref. [1] because (i) here B 2 is is used and (ii) a relative rather than absolute pressure 
is used. The definition in Ref. [1| used, in contrast, the average poloidal held linking the 
circular surface lying in the z = 0 plane between the geometric axis and the magnetic axis. 
Because the definition of /3 re i uses the relative hydrodynamic pressure, it is seen that /3 re i can 
be positive or negative. In particular, if P ax i S is smaller then P^ c , then f3 re i will be negative. 
The definition of /3 re i is useful because it provides a simple mathematical way to distinguish 
toroidal equilibria with inverted (3 profiles from those with normal (3 profiles. The former 
are toroidal magnetic bubbles while the latter are toroidal confinement configurations such 
as spheromaks and tokamaks. 

On expressing the magnetic held as 


B = — (V0(r, z) x V</> + z)V</>) 

Z7T 


( 6 ) 


where / = ‘InrB^j is the poloidal current, MHD equilibrium JxB = VP can be ex¬ 
pressed as the Grad-Shafranov equation AB 


d f 1 dtl> 
dr \r dr 


d 2 i\) 9 9 dP 2 dl 

_ + 4x/Ior _ + Mo/ _ = 0 . 


( 7 ) 


We assume that P is a linear function of the poloidal hux 0 and so can be expressed as 


P = 


Pn 


Ic 


■i> 


Paxis^lc Plc^ a 


( 8 ) 


('Tpaxis fplc) (jPaxis ” 0 Ic ) 

where ipi c is the last closed hux surface of the configuration. 

The poloidal current is similarly assumed to be a linear function of the poloidal hux and 
can be expressed as 

Po I = A0. (9) 
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We note that the assumed linear dependence in Eqj9] differs from the assumption used in 


Solov’ev-type solutions such as in Ref. [4] where it is assumed that I 2 ~ ip + const. For 
the linear dependence assumed here, Idl/dip is linear in 0 whereas for the Solov’ev-type 
assumption, Idl/d 0 is a constant. 

Using EqlUl the toroidal component of EqJH] gives 


B, 


\i> a 


2nr n 


( 10 ) 


The gradient of P with respect to 0 can then be expressed in terms of (3 re i as 


d2 

^ axis 


d_P = ^ _ 

dip dop^axis (1 Iplcj^Paxis) 

Since B po i vanishes at the magnetic axis, B 2 xis = B 2 axis and so 

dP _ g A Ipaxis 

dip rel 4n 2 rl xis no (1 - ipic/ipaxis )' 


( 11 ) 


( 12 ) 


III. CYLINDRICAL SOLUTIONS TO FINITE /3 rei GRAD-SHAFRANOV EQUA¬ 
TION 


We now introduce dimensionless quantities 


ip 


’axis 


r 



T axis T* axis 


A r 


axis 


SO EqJT] can be expressed as 


We define 


_d_ 

dr \r dr J 


d 2 ip 
dz 2 


+ A 2 



o. 


PrelT 
(1 - tic) 


(13) 


(14) 

(15) 


so Eq jm becomes 


_ y /i«\3 + X2 

dr \r dr ) dz 2 


0. 


(16) 


We assume a solution of the form 


X = rg(r) cos (kz) 


(17) 
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so Eg JI61 becomes 


d 2 g 1 dg (_ 2 1 \ 

ft5 + FaF + E ~^) g -° 

where 

f = A 2 - k 2 . 

Equation [IS] is Bessel’s equation with general solution for real 7 


(18) 

(19) 


g(r) = + oy\\{-ir) ( 20 ) 


where ctj and ay are constant coefficients to be determined by boundary conditions. 
From EqsflSl and [171 the solution to the normalized Grad-Shafranov equation is 


ip = r (ajJipyr) + ayY^r)) cos (kz) — r 2 Q 

(21) 

where 

^5 drel 

^ ~ l-'lplc 

(22) 

However, ip — 1 is required at r = 1, z — 0 (i.e., at the magnetic axis) so 

<JjJ\{i) + <JyY x {i) = 1 + Q . 

(23) 

The following three Bessel identities where C n = J n or Y n will now 

the rest of the discussion: 

be used repeatedly in 

dc f s) = ft M 

as 

(24a) 

2 (scyo)) = s c„( s ) 

(24b) 

sG 2 (s) = 2C\(s) — sCq(s) . 

(24c) 

The magnetic axis is also where dip/dr vanishes and so taking the derivative of Eql2T1 with 

respect to f. usinef Eql24b and then setting f — 1 and z = 0 elves 

2 

ffi/.M7) + 0V*o(7) = -Q- 

7 

(25) 


Equations [23] and [25] constitute two linear inhomogeneous algebraic equations for the coef¬ 
ficients crj and cry. Solving these equations for aj and ay and using the Wronskian 


Ji(s)Y 0 {s) - J 0 (s) YM) = — (26) 
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and Eq J24cl gives 


u = y (>0(7) - QY2 (7)) 
<rr = y (~Mi) + QMl)) 


(27a) 

(27b) 


IV. SPHEROMAK-TYPE SOLUTIONS 


Spheromaks are singly-connected Grad-Shafranov equilibria (i.e., there is no “hole” in 
the “doughnut”) and so the domain includes f = 0. A spheromak therefore cannot contain 
a Y\ (yf) component because Yi(tt) diverges at f = 0. It is thus necessary to impose cry = 0 
for a spheromak in which case Eq l27bl yields the relation 

Ml) 


Q = 


Mi)' 


Substituting for Q in Eq l27al and using Eqs l24cl and [26] gives 

2 


0J = 


(28) 


(29) 


7^2(7)' 

Using Eq l24cl to substitute for ^( 7 ) in Eq 1251 shows that EqJ25]can alternately be written 


as 


so one can also write oj as 


1 M 1 ) 

l + Q 
■A (7) ‘ 


oj = 


Because cry = 0 for a spheromak Eqs[22] and [28] show that a spheromak has 

A* = (1 - M) 

Mi) 


and 


<?j = 


1 + 


firel 


(30) 

(31) 

(32) 

(33) 


Ji(7) V 1 -iplc 

On substituting for <jj and Q in Eqj21] the solution to the normalized Grad-Shafranov 
equation becomes 


ip = r 


-Jli.l'f) (\ fir el \ /J _2 (dr el 


1 + 


cos (kz) — r 2 


(34) 


Mi) V ' 1 v ’ 1 -i’k 

If /3 re i = 0 and i/Sj c = 0 are additionally assumed, the standard result for a zero-beta 
spheromak in a cylindrical flux conserver of radius a is retrieved, namely 7 = x 0 i = 2.405 
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where xoi is the first root of Jq. Since 7 = 7 r ax is, and the last closed flux surface is at the 
cylinder radius, then the assumption fii c = 0 and /3 re i = 0 in Eq. [3T1 implies ^ 1 ( 70 ) = 0 in 
which case 7 a = x u = 3.83 where in is the first root of J\ . Thus, for a fi re i = 0 spheromak, 
Taxis/& = 7 /( 70 ) = x m /xu = 0.63 as is well known. Equation 1371 shows that spheromaks 
with finite positive /3 re i are restricted to the range 0 < 7 < 2.405 but, as will be discussed 
in Sec lYll physically relevant non-spheromak configurations with negative /3 re i exist when 
7 > 2.405. 

Substitution of Eqj32] into Eqj34] gives 

1 / 2 f - \ 

^ = Mfi) \J Jim cos ( k *) - (35) 

which reverts to the fi re i = 0 solution when 7 = 2.405 as can be seen using Eq l24cl to give 
7 ^ 2 ( 7 ) = 2Ji( 7) if J 0 ( 7 ) = 0. 


V. SAFETY FACTOR OF SPHEROMAKS WITH FINITE fi rel 


The last closed flux surface of a spheromak has i/>i c = 0 and Pi c = 0 in which case Eqs. [5] 
and m give 

(36) 


firei = 4vr r axi 


2 2 h'O Paxis 


axis A 2 ?/; 2 

/ * r axis 


and Eq 1341 becomes 


V 1 = r ~T7~V ( 1 + Prd) COS ( kz ) “ PrelT 2 


(37) 


which is the same as Eq.(2) of Ref. lj except for the different definition of fi re i. 

In order to determine q ax is , Eqj3] shows that it is necessary to calculate /do J<f>,axis/ 
Equation [2] shows that 


V’a 


j T Uj Ju ta 

lk) ' * = ~ 2 vrfr 3 

1 'axis 

so, using Eq JlOland EqR4l it is seen that 

ho '7,6 ,axis 


_ d f 1 di/j \ d 2 ^ 

df \ f dr ) dz 2 


— (1 + (3 re i) A. 


(38) 


B(f>,axis 

Thus Ho Jfraxis/B^axis = A only if (3 rel = 0. Inserting EqSSlin Eq|3] gives 

e 1 ' 2 + e-V 2 


(39) 


Qaxis 


A (1 + fir el) 


(40) 
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f.fl 


which differs from Eq.(30) of Ref. [l] by having an extra and important factor of (1 + f3 re i) 
in the denominator. 

From Eqm and use of the Bessel identities it is seen that 

$ zz ) axis = 

(M) axis = 

so the cllipticity is 

e = 

This indicates that the poloidal flux surfaces will be circular near the magnetic axis (i.e., 
have e = 1) if 7 = k in which case A = s/2'y. Combination of Eqsl32l 1401 and 02] gives 


2 k 2 J x ( 7 ) 

7 Mi) 

(41a) 

-2iMi) 

Mi) 

(41b) 

7 2 

k 2 ' 

(42) 


Qaxis 


A 


A J 2 ( 7 ) 


(43) 


^kl + /3 re i 2k Ml) 

Equation 021 has been validated by direct numerical integration of field lines in the vicinity of 
the magnetic axis of a magnetic configuration characterized by Eqj 6 ] with -0 given by Eql35l 
In the (3 re i = 0 limit, J 0 ( 7 ) = 0 and q axis —* A/ (yfc) which is Eq.(33) of Ref.[ 1 ], but for finite 
positive /3 re i, Eq. 03] shows that q n ^ is reduced from its /3 re i = 0 value. 

The safety factor at the wall is 1] 

Qwaii = TTgr cos_1 (Ml)) ( 44 ) 

and so the ratio of safety factor at the wall to that at the axis is 

qwaii Ji(i) cos" 1 (J 0 (7)) 

— = - ( - 45 ) 

Qaxis J 2 \)) 7T 

which is plotted in Fig.l. Contrary to Ref.[1] it is seen that the shear (difference between 
q wa ii and q ax is) decreases with increasing f3 re i (i.e., with 7 decreasing below 2.405). Using 
16 = In (cos 6 + i sin 6 ) to write 


cos 1 (J 0 (7)) = —* 111 (-Ml) + “ (-^oW ) 2 


(46) 


and then using Jo(t) = 1 — 7 2 /4 for 7 <C 1, it is seen that for 7 <C 1 

r ' 2 . i \ 7 

1 - I - 

4 ' V 2 J ~ V 2 ' 


cos 1 (Jo ( 7 )) — — i hi ( 1 — + i —p ) — —j=. 


(47) 
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(48) 


Since ^ 1 ( 7 ) — 7/2 and 7 ) — 7 2 /8 for 7 -C 1 , Eqj45]has the limiting behavior 

— -a —t- = 0.900 for 7 < 1 

Qaxis 

which is seen in Fig. 1. Furthermore, Eqj32]has the limiting behavior 


Prei ->• vy for 7 -C 1 ; 
7 

i.e., /3 re / diverges at small 7 which is also seen in Fig. 1. 


(49) 


We note that numerical calculations reported in Ref. 


8 ] assumed / 2 ~ 0 2 (1 + 2a0/3) and 


dP/dil> ~ 0 — 0o hi a spherical geometry and found that the gradient of the shear had a 
strong dependence on a. The analytic solution given here would correspond approximately to 
the a = 0 numerical solution reported in Ref.[ 8 J; the correspondence is not exact because of 
the different assumptions for the dependence of P on 0, the shape of the boundary (cylinder 
v. sphere), and the assumption of a central hole in Ref. 8 ]. 


VI. TOROIDAL MAGNETIC BUBBLE: NEGATIVE /3 re i 

We now consider the situation where (3 re i < 0 and 0/ c 7 ^ 0. We consider the ay = 0 case 
first as was assumed for spheromaks and then later consider the more general case where 
both <jj and ay are finite. 


A. oy = 0 case 

In the oy = 0 case 0(f, z) is mathematically identical to the spheromak solution consid¬ 
ered in Sec.IV, i.e., Eq [35] provides the relevant flux function. The difference here is that 
iJji c is no longer assumed to be zero. Plots of i/>(f, z) using 7 > 2.405 show that ^(f, 0) has 
periodic maxima and minima because of its -J\ ( 7 f) dependence. Equation [13] defined -0 to 
be unity on the magnetic axis, i.e., 0 (f, 0 ) = 1 at f = 1 and the magnetic axis was defined 
to be where 0 was a maximum or minimum. Because of the oscillatory behavior of Bessel 
functions, maxima or minima of 0 occur not only at r = 1 but also for f > 1. However, 
the maxima and minima occurring where f > 1 do not have 0 = 1 and so do not satisfy 
the 0=1 condition given in EqIT3l Thus, only the maximum of 0(r, 0 ) at f = 1 will 
be considered since maxima or minima at larger f do not satisfy the 0=1 requirement 
stipulated in EqJl3l 
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Examination of Eqj35] shows that $ is independent of z if = 0; at this radius 

f = xn/'y where X\\ = 3.832 is the Erst root of J\. We now show that this radius f = Xu/'J 
is infinitesimally larger than the radius of the last closed flux surface. Since xjj is independent 
of z when = 0 , the flux surface passing through f = xn/ 7 , z — 0 must be a straight 

vertical line, i.e., ^(xu/^f,z) = ^(xu/ 7 , 0 ) for all z. Because a straight vertical line goes to 
z = ± 00 , the flux surface passing through f = xn/ 7 , z = 0 is open. Immediately to the left 
of this line the flux surfaces are closed and so the last closed flux surface is at the radius f; c 
where 


ric = lim 
< 5->0 



Xu 

7 ' 


(50) 


This can also be seen graphically from the flux surface contours shown in Fig.2 (to be 
discussed in more detail later) where it is seen that a straight vertical line separatrix lies 
between the blue-purple closed flux surfaces having magnetic axis at f = 1 , z = 0 and the 
green-orange flux surfaces to the right. Equation [501 gives the radial location of this vertical 
line. 

A toroidal inverse aspect ratio (ratio of torus minor to major radius) can be defined as 


Tic T axis — 

£ = - = Tic 


1 = 


xu - 7 
7 


(51) 


Using Jii^jric) = 0 at the last closed flux surface, Eq. I35lmav be evaluated at r = ri c , z = 0 
to give 


4>ic = 


%u Mi) 

7 2 -Mt)' 


Inserting in Eql32l gives 




(52) 


(53) 


7 2 ^ 2 ( 7)7 Ml)’ 

In order to have r; c > r axis Eql5H shows that it is necessary to have 7 < Xu = 3.832. A plot of 
Eq J53l shows that f3 re i is negative if 2.405 < 7 < 3.736; /3 re i changes sign at 7 = 3.736 because 
the quantity in parenthesis in Eql53l changes sign at 7 = 3.736. Thus if 2.405 < 7 < 3.736, 
firei is negative and also r ic > r axis . 

Because the minimum of (1 + s)s occurs when s = —1/2, identifying s = 

x ii^o( 7 )/ ( 7 2 ^ 2 ( 7 )) it is seen that (3 re i is at a minimum when x\ x ^o( 7 )/( 7 2 ^( 7 )) = - 1/2 in 
which case 


min [j3 re i] = 


7 


4x? 


(54) 
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Using the Bessel identities, the magnetic field components are 


B r = 


Praxis d'lp Tpaxis 2/c J\ i^T ) 

2 ^Lm 7 -Mi) 


Ba, — 


B y = 


2 *r 2 axis r dz 

A”0 'Paxi: 


sin(A:z) 


2tt r 

Ipac 


2nrt 


dpi ip a 


Mi) 


2nr 2 ,s dr 2nri 


—Mir) cos (kz) — r J 0 ( 7 ) 
7 


— ( Ml r ) cos (kz) - Ml)) 


Ml) 


Using Eq J24cl and Eq 15 5b I it is seen that 

B(j)^axis 




2n r 2 


A. 


(55a) 

(55b) 

(55c) 

(56) 


A normalized magnetic field can be defined as B = B(r, z)/B ( j ) ^ axis with components 

(57a) 
(57b) 


2 k Mir) . T 

B r (r,z) = — j sin (fez) 


7 A Ml) 

^Mir) cos(kz) — fMl) 
BMi*) = ’ 


B z (r,z) = 


Mi)_ 

2 (Mir) cos(kz) - Mi)) 


(57c) 


A Ml) 

As required, both B r and B z vanish on the magnetic axis (i.e., at f — 1 , z — 0 ) and B $ = 1 
on the magnetic axis. 

Equation [351 with 2.405 < 7 < 3.736 thus gives the flux surface for a magnetic bubble, 
i.e., a toroidal configuration with closed field lines where the pressure on the magnetic 
axis is lower than the pressure at the surface of the toroid. The direction of the J x B 
force is thus outwards rather than inwards in contrast to a tokamak. This configuration is 
relevant to axisymmetric interplanetary magnetic clouds ejected from the sun by coronal 
mass ejections. Spacecraft measurements indicate that P is smaller in the interior of these 
clouds than outside so these clouds have negative Mi- Another possible situation would 
be in the solar interior where a toroidal bubble configuration as described here would be a 
toroidal region of stronger magnetic field but reduced hydrodynamic pressure compared to 
the surroundings. 

As a concrete example of such a configuration, consider the situation where 7 — k — 2.5 
and A = \f2l- In this case e = 1 so the poloidal flux surfaces are circular near the magnetic 
axis, the last closed flux surface is at M = 0.25 and from Eqj53]/3 re ; = —0.081. From EqlSTl 
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it is seen that the inverse aspect ratio is e = 0.53. Figure 2 plots contours of ij){f , z ) and it 
is seen that the last closed flux surface intersects z = 0 to the right of the magnetic axis 
at indeed ?y c = Xn/'y = 1.53. Figures 3, 4, 5, and 6 plot ^(f, 0), B^f, 0), B z (f, 0), and 
B 2 (r, 0) respectively. 

From Eqjn] it is seen that 

/fo p _ p 1 /a 

7^2 ^ arris 7^2 lc “r l^rel 

^axis ^axis 

so the hydrodynamic pressure on the magnetic axis is lower than on the last closed flux 
surface. If P ax is is set to zero, then the external pressure would be 

= -Prel (59) 

in which case the configuration would be a vacuum at the magnetic axis (zero plasma 
pressure) with increasing pressure going away from the magnetic axis toward last closed flux 
surface. 

If 7 is further increased, the sign of (3 re i can become positive again in which case the equi¬ 
librium will become tokamak-like (higher pressure on magnetic axis). Additional increase 
of 7 will cause (3 re i to oscillate in sign giving a sequence of bubble-like and tokamak-like 
configurations. Also, for a given configuration one could elect to truncate the flux at some 
value larger than ipi c and so obtain a smaller aspect ratio equilibrium. In accordance with 
the Shafranov virial theorem, any one of these configurations will involve a jump in the 
magnetic field at the surface of the toroid if it is assumed that at the surface the external 
magnetic field differs from the internal field. This jump corresponds to the existence of 
surface currents. I 11 a tokamak these surface currents are provided by a set of coils im¬ 
mediately external to the toroidal volume and these coils are called the vertical held coils. 
The held produced by these coils is mainly in the z direction and will be referred to here 
as B e z xt . This held B ext constitutes a portion of the total held inside the toroidal volume 
and provides equilibrium in the major radius direction. This takes place via a radial force 
~ J,pB e z xt directed towards f = 0 that balances the radially outward hoop force as well as 
some hydrodynamic pressure forces. The hoop force is a property of any toroidal current 
system and occurs because a toroidal current produces a stronger poloidal held near f = 0, 
z = 0 (inside) than at f = f/ c , z = 0 (outside). This stronger poloidal held on the inside 
compared to the outside corresponds to greater magnetic pressure on the inside than on 
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the outside; for low f3 the force resulting from magnetic pressure imbalance dominates any 
hydrodynamic pressure imbalance. Without the offsetting force provided by B z xt the hoop 
force would act to expand the torus major radius. 

At first sight it might appear that the flux contours in Fig .2 are such that the magnetic 
pressure is higher on the outside than on the inside because the midplane poloidal flux 
surfaces in Fig. 2 are more tightly packed outside the magnetic axis (e.g., at f ~ 1.5) than 
inside the magnetic axis (e.g., at f ~ 0.25). However, the density of held lines and hence 
the poloidal held is nevertheless stronger inside the magnetic axis than outside because of 
the inverse f dependence in B z = (2nr)~ l d'ili/dr. The twice as tight midplane hux surface 
packing in Fig. 2 at f ~ 1.5 compared to at f ~ 0.25 gives a twice as large d^/dr on 
the outside compared to the inside. However, this twice as tight radial packing is overcome 
by the (27rr)" 1 factor, a toroidal geometry effect that produces an approximately six-fold 
inside-to-outside enhancement with the net result that \B Z \ is about three times larger at 
f ~ 0.25 than at f ~ 1.5. This three-fold inside-to-outside ratio of \B Z \ is evident in Fig. 
5. 

In order to have the B z xt required for equilibrium, it would be necessary to have surface 
currents flowing on the surface of the toroid. Since there are no powered coils to sustain 
surface currents exterior to a magnetic cloud, it is unlikely that such surface currents would 
exist in the magnetic cloud context. Without the B e z xt provided by surface currents (and 
intrinsic to the equilibrium given here), the hoop force resulting from the imbalance between 
B 2 on the inside and B z on the outside will cause the major radius of magnetic clouds to 
increase with time. The difference between poloidal flux surfaces with and without incorpo¬ 
ration of B z xt is of the order of the inverse aspect ratio £ because B z xt is a toroidal effect 
and so scales as e. 

B. Finite crj and ay case 

The spheromak solution required cry to be zero to avoid singularity at f = 0. The magnetic 
bubble solution discussed above used the same functional form as the spheromak solution 
(i.e., had cry = 0 and used Eql34l) and found that a tokamak-like solution with /3 re i < 0 (i.e., 
inverted beta profile) occurred if 2.405 < 7 < 3.736. If f = 0 is excluded from the domain so 
the configuration is doubly-connected, the singular nature of Y\ (yf) at f = 0 is no longer a 
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constraint and the more general solution given by Eqs. EU I27al and I27bl can be used. The 
consequence of imposing cry = 0 was for Eq J27bl to force the relationship between Q and 7 
given by Eg EM If oy is not forced to be zero, then this relationship between Q and 7 is no 
longer imposed and the only remaining condition is that the domain must exclude f = 0 . 

Consideration of EqEUand recalling the discussion that led to EqEtUI shows that ^(r^, z) 
is independent of ! at fi c where fi c is now defined by 


vjJi(iri c ) + ayY^ric) = 0 . 


(60) 


Thus Eq 160] provides a radial shift of the location of the last closed flux surface and generalizes 
the discussion that led to Eqj50] Because <jj and cry depend on 7 and on Q (hence on /3 re i ), 
Eql60l shows that fi c depends on both 7 and (3 re i . However, by assumption f; c > 1 (last 
closed flux surface radius is to the right of the magnetic axis) which restricts the allowed 
values of 7 and /3 re i. Introduction of the Y\ (jf) solution and the coefficients 07 and cry is 
thus analogous to generalizing the solution of some harmonic equation from being sin (kx) 
to being sin(fcx + A) = sin(kx) cos(A;A) + cos(kx) sin(fcA) where sin (kx), cos (kx) are the 
analogs of (qr), (qf) and 07 , 07 - are the analogs of cos(/cA), sin(/cA). Introducing hnite 
A changes the phase of the solution and shifts the location of the solution. 

Substitution for 07 and ay in Eqj60] using Eqs J27al I27bl gives 

Y 0 (7)Ji(Wic) - Jo(7) Y i(Wic) 


Q(7, Ac) = 


( 61 ) 


Y 2 ( 7 ) Ji(7Ac) - MrfYii'jric)' 

It is seen that Eq J6T1 reduces to Eqj28]if J\{jfi c ) = 0, i.e., the situation considered in Sec JVI Al 
and that Q becomes infinite when ?y c is such that the denominator in the right hand side 
of Eqj 6 T] vanishes. 

Using Eqsj60]in Eqj2]]it is seen that the last closed flux surface is given by 


V’fc = ~rlQ 


(62) 


and inserting this in Eql22l gives 

(3rel — Q + rfcQ 2 ■ (63) 

The derivative of Eq]63l shows that the minimum possible (3 re i is f3 re i = —1/(4 rf c ) which 
occurs when Q = — l/(2ff c ); this generalizes Eq]54l 

Figure 7 plots the dependence of Q and (3 re i on f; c for 1 < U c < 2 with 7 = 2.5; it 
is seen that, as predicted, (3 re i has a minimum at /3 re i = — l/(4f/ c ) which occurs when 
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q = -i/(2f,y. It is also seen from this figure that when n c = £11/7 = 1.5328 the Sec IVI Al 
result /3 re i = —0.081 and Q = 7) = —0.108 is recovered. For this 7 = 2.5 value 

the denominator in Eq ffTTI vanishes when f/ c —> 1.885. 

The following chain of dependence thus exists for doubly-connected configurations: 


1. Independent values for ri c and 7 can be selected which then determine Q via Eq. [6TJ 

2. Using Eqs l27al and I27bl in EqEUthe flux function 0 ( r,z ) is given by 


//- -n 7T--J Oo(T) -<3^2(7)] . / r-\ _ 2n 

0( r , zj = — 7 r < > cos(kz) — r Q, 

2 ' + [-^ 0 ( 7 ) + QMi)l Yi(V) 


(64) 


3. /3 re i is given by Eql63l 

4. {j)i c is given by Eql62l 


This chain of dependence for doubly-connected configurations differs from that of a finite 
P spheromak. Specifically the chain of dependence for a finite /3 spheromak is: = 0 is 

imposed because of the singly-connected topology, f/ c is determined from setting the left 
hand side of Eqj35]to zero on the midplane, and EqJ32l gives (3 re i = '/oOtV'M'?)- 

Another and equivalent point of view differentiating singly- and doubly-connected con¬ 
figurations from each other is the following: 

(i) because the midplane of a singly-connected configuration contains f = 0 and because 
"0 = 0 at f = 0, the last closed flux surface for a singly-connected configuration must always 
have 0 = 0, 

whereas in contrast, 

(ii) because the midplane of a doubly-connected configuration excludes f = 0, the last 
closed flux surface of a doubly-connected configuration cannot be 0 = 0 as such a flux 
surface would have to pass through f = 0. 

The magnetic field components associated with Eql64] normalized to B 0 /lxis = 
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^•4^ axis / (2 ^Lis) are 


B r (f, z ) 
z) 

B z (f,z ) 



[>0(7) - < 5^2 (7)] -Mtt) 

+ [-Jo (7) + <5^2(7)] Y^f 

[>0(7) - < 5^2 (7)] Mir) 

+ [-Jo (7) + <3^2(7)] *1(7*0 

[>0(7) - QY2 (7)] Mir) 

+ [~Mi) + QMi)) Mir) 


sin(fc2) 


cos (kz) — rQ 


cos (kz) — 2 Q 


(65a) 

(65b) 

(65c) 


Using Eqs J24cl and 1261 it is seen that Eq 1651 reverts to EqJ57l when Q = Jo(i)/Mi)- 
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FIG. 1: q wa ii/qaxis plotted as black solid line v. 7 from Ecil45l and (3 re i plotted as red dashed line. 



0.0 0.5 1.0 1.5 2.0 

FIG. 2: Contour plot of V>( r,z ) given by Eq|35]for 7 = k = 2.5. The radius of the last closed flux 
surface is at f = 1.53. The hydrodynamic pressure is lower in the pink region than in the green 
region so the configuration is a magnetic bubble. 
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FIG. 3: Plot of 0) as given by Eq[35] for 7 = k = 2.5 and A = 2 ^/ 7 . The radius of the last 
closed flux surface is at f = 1.53 and t/>j c = 0.25. 



FIG. 4: BAr, 0); Note that magnetic axis is at r = 1.0 and that maximum of B^ occurs to left of 
magnetic axis. 



FIG. 5: B z (r, 0); note that B z changes sign at magnetic axis. 
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FIG. 6 : Bj. + + B% as function of r for z = 0. 



FIG. 7: Figure 7. Blacksolid line is plot of Q( 7 , r; c ) as function of r/ c with 7 = 2.5; red dashed line 
is plot of Prei- For this 7 , Q diverges when f/ c = 1.885. 


21 








